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Abstract 

The role of lattice asymmetry parameter £ in description of the SU(2) gluodynamics 
phase structure at finite temperature is studied analytically. The fact that renormal- 
ization group relations which permit to remove the lattice asymmetry parameter from 
the thermodynamical quantities in the "naive" limit don't do the same in the approx- 
imation SU(N) ~ Z(N) keeping the "naive" limit results at the same time was point 
out. An additional condition which would fix £ is needed for this dependence removal. 

INTRODUCTION. 



Following the heart point of the renormalization group theory after any 
regularization removal (in particular, lattice one) observed quantities should 
go to their physical values and should not depend on the renormalization 
scheme. So, although going to the continuum from a symmetrical lattice 
and from asymmetrical one are two different renormalization procedures in 
the gauge fields theory it is natural to expect that in the continuum limit 
no lattice asymmetry effects will be revealed. At the same time, strictly 
speaking, the renormalization group trajectories are universal (i.e., they don't 
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depend on the renormalization scheme) only near the critical point at zero 
coupling while at intermediate and strong couplings this statement is still 
questionable. It is quite probable the trajectories will depend on the degree of 
lattice asymmetry Moreover, today's level of this problem understanding 
does not permit to exclude other possible revealings of the lattice asymmetry. 

In this paper we consider the finite temperature SU (A^)-gluodynamics on 
an asymmetrical lattice 

a^^a v \ /i = 0, 1,2, 3; a N = [3; a k N k = L k ; k = 1, 2, 3. 
a M - the lattice spacing along the direction fi. The most wide-spread ex- 
ample of such asymmetry introducing is the Hamiltonian limit |T], 0] (where 
the lattice distance along the space directions is fixed and the spacing along 
the time direction is going to zero). The renormalization procedure inde- 
pendence actually means that renormalized couplings of both theories (for 
example, Hamiltonian and Euclidean SU (N) lattice gauge theories) should 
be get equal. This naturally imposes the relation on the scales j^. Following 
chosen approach [[I], |2j] this is enough for two renormalized theories to give 
the "same physics". Indeed, if in the standard Wilson's action 
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X,fl<V N c 



S = /V £ ReTr(l-— U ll {x)U v {x + n)U+{x + v)U+{x)) {V 



we restrict ourselves to only "naive" limit (which assumes the expansion of 
SU (N) matrices around I-matrix, close to I-matrix or any fixed matrix) 

U,(x) = e ia ^ A ^) ~ i + ia^A^x) (2) 

and we take 

_ 2N aoa^as . . 

K>(j,v — 2 2 v ) 

then at least in the Hamiltonian (H) and the Euclidean (E) [a v = a; g% —> g 2 E ] 
limits the spectra of hamiltonians coincide. Moreover, not only in the classical 
case (g 2 E = gjj), but in the quantum case (g E ~ g\ + O(gjj)) too. 

However, investigation of the finite temperature SU (N) lattice gauge the- 
ory in "naive" limit is carried out using SU(N) matrices expansion around 
fixed matrix thereby leaving the center Z(N) out of the game. But the role of 
the gauge group center in the finite temperature QCD phase structure is well- 
known. Therefore, we will study the lattice asymmetry effects contributed 
by the center Z(N). To do this we will use the approach SU(N) ~ Z(N), 
i.e., we will consider the finite temperature lattice gauge theory in which the 



2 



link variables are the center elements of the SU (N) gauge group 

a fl (x) = e—; ^ = 0,1,2. ..N-l (4) 

On the other hand, for theories with gauge groups which have a trivial 
center (for example, SU(N)/Z(N)) it may be quite reasonable to restrict 
ourselves to "naive" limit, whereas for the SU(N) gauge theories such a 
restriction requires special substantiation because there is no naive limit for 
the center Z(N) in common sense. Building the continuum quantum field 
theory becomes possible only if Z(N) lattice gauge system undergoes the 
second order phase transition at zero of the renormalization group function 

b a a. 

To ensure agreement with the "naive" limit results we chose (|3|) to 
guarantee the coincidence of naive limits for the models with SU (N) gauge 
group on the symmetrical lattice and asymmetrical one. 

In the future we will calculate the corrections to our approximation. As 
a justification of our approach we should note that effect of the quantum 
fluctuations, which are taken into account besides the Z(N) topological ex- 
citations, amounts to simple change in the coupling constant [@, 0, |7f 

Knew = K i d - {N 2 - l)/4 (5) 

Additional term (|5]) depends neither on a M nor on k,^. Therefore, it does 
not depend on the lattice asymmetry. So, there is a reason to hope that 
including the quantum fluctuations besides of Z(N) solutions in our theory 
on asymmetrical lattice will not change ([5]), i.e. will just result in the finite 
renormalization of couplings. 

It will be enough for our purposes to consider the simplest case of lattice 
asymmetry, when k,q h — k t and K nm — K a , i.e., we allow the couplings to be 
different for plaquettes lying in the temporal and in the spatial planes. Our 
results may be easily generalized to the case (§). We will also limit ourselves 
to the analyses of four limiting cases: 

1) K a — > 0; 2) k t — > 0; 3) K a — > oo; 4) k t — > oo; (6) 

(the estimations of phase structure in the whole plane K T 0K a we are planning 
to present in next paper). 

As known our Wilson's action is splitting into two parts: 

- S(a) = n a - -^Re Tra n (x)a m (x + n)cr+(f + m)cr+(f)^ + 
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«t 2^1 - ]y Re Trcr °(^' T ) a n(x, r + 0)(Jq(x + n, r)cr+(f , r)J (7) 

where k ct = 2N a T j g\a a and ft T = 2N a a j g 2 T a T \ a a (a T ) is the spatial (tempo- 
ral) spacing and the sum P a (P T ) runs over all purely space-like (time- like) 
plaquettes. 

Let's first consider j a = th.K, a and 7 T = th.K T independent temporary 
ignoring the dependence of g a (g T ) on a a (a T ) through the renormalization 
group relations, i.e., supposing that by varying a a and a T we can get to any 
point of the square (0 < 7 r < 1; < 7 CT < 1). We will investigate each 
of the limiting cases separately, joining the results for the couplings critical 
values through the same parameter. 



THE AREA OF SMALL n a . 

In the case 7 ff ~ we may throw off the magnetic part of action thereby 
leaving only the spin configurations with time-like plaquettes. In the static 
gauge we can sum over the space spin configurations {cr n }, finally gaining the 
Ising model which is much known of. 

Now let us dwell upon that. Fixing the static gauge o"o(x, r) = u$ we 
consider the following gauge transitions 

<r+(x , r) -> u^ n a+(x, t)uI (8) 

Imposing the periodical boundary conditions in temporal direction upon 
cr-matrices: a n (x, t = 1) = a n (x,r = N T ) , all matrices u$ are grouped on 
the last links resulting in the Polyakov's loops 

an = n <o* = ( 9 ) 

T=l 



We come to the following action 



N T -2 

E (l-iReTia n (x,r)a+(x,r+l)) + 



r=l 



(1 - ^Re TvQ,a n (x, N T - l)n± +n a+(x , 1)) ] (10) 
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Summing over all configurations {cr n } in Z(2) gauge system we get 

- s Pr = K \j:n^ +n (ii) 

x,n 

th« T = (thK T ) Nr = 7r 

As known, for the Ising model there exists the critical value of the coupling 
<-j7 T = 7 C separating two phases. Consequently, the partition function fllPl) will 
have a critical point at 7^: = jI' Nt - 



THE AREA OF SMALL k t . 

In the other limiting case 7 r ~ the functional integral is highly peaked 
about configurations with space-like plaquettes and we may throw off the 
electric part of action. In this case the partition function turns out to fall 
to N T equivalent disconnected contributions each of which is from a separate 
3-dimensional layer. We get the set of standard 3-dimensional Wegner models 

Sp a -EE /V<7 n <7 m <r;<T* (12) 

T= l x,n,m 

There is no interaction between the layers, so summing over o~ n (x : r) may be 
done independently for every r = const layer. 

The lattice dual to original hypercubical lattice can be constructed by 
shifting the lattice by half a lattice spacing in each direction (see, for example, 
H). Geometrical duality transforms g-dimensional manifolds into (d — tri- 
dimensional ones. "Dual coupling constant" R for the Z{2) gauge theory 
K a = — |lnth« CT is a monotonically decreasing function of the "original" 
coupling. The set of 3-dimensional Wegner models transforms into the set of 
3-dimensional Ising models with spins in sites under duality transformations. 
These Ising models exhibit the transition at 7 C simultaneously. Therefore, for 
the partition function ( |I2|) there is a critical point at 7^ = jt^- 



TWO ANOTHER LIMITING CASES: n T -> 00 AND n a 00. 
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When considering the duality transformation in 4-dimensional space-time 
it should be pointed out that just space- like plaquettes transform into time- 
like ones. In the other words, 

< = -^lnth« T or = (13) 

2 1 + 7r 

and vice versa, k! t = — ^lnth«; CT . 

This statement becomes clear from the following. Let us rewrite the par- 
tition function of Z{2) system in the form 

g = q k ij.v Y^x,tiv u x,ixv (14) 
M 

= E II ch/^(l + Dx^thK^) 

{a} 

= (""En E(a^th« M „)^ 

= e-^E n(e lnthK -)^ II E (a)^-^^ 

1^,1. Zinfcscr=±l 

3 



e 



-TV/ 



e n (e^y-^ 11 n 2^ 2 ( e ^ ±1 ) 



where U x ^ v = <7 li (x)a 1/ (x + ^)cr*(x + v)al(x) and 
f = NrNlT.^thK^. 

We introduced a new set of variables {q} - one for each plaquette. The 
partition function is not equal to zero only if q x ^ v satisfies the following 
condition on the sum over six q x ^ v (associated with six plaquettes which 
adjoin the link x, v, see fig. 1) : 

I 3 3 

« E {Qx,iw + 1) = mo d2 or E £r,/^ = 2 mod4 (15) 

The solution of last equation can be found if we associate every q XjjJLV with 
one of the cube plane and 

qx,nv = s p (x)s\(x + p)s*(x + X)s* x (x) 

v ± p^p^X (16) 

where the dual link variable s p {x) is the element of the Z(2) group. It becomes 
intuitively evident, if we consider the starting case when all s are equal 1. 
It dictates for T,^±v Qx,fj,v to be equal 6 mo d4 = 2 mo d4- Every link enters the 



solution twice (because the plaquettes form a cube) and changing the sign of 
a link to opposite results in changing £^ = _ 3 . M ^ ±l , q x ^ v only by ±4. 
Consequently, in the plane 

Kit ® Kfj there is the self-duality line (see fig. 2). 
R. Balian, J.M. Drouffe, C. Itzykson [J9j pointed out the possibility of the 
critical behaviour at k c = 0.44 for the 4-dimensional Z(2) pure gauge theory 
on symmetrical lattice supposing this critical point is single. 

So, now we realize that under duality transformation our original 4- 
dimensional theory (0) transforms into the same one but with new coupling 
constants. 

-S' = + (17) 

For this dual representation we can consider the two limiting cases in 
precisely the same way we have done for the original theory. 

The case j' T = th.K,' T ~ for the original theory means 7 CT ~ 1 in according 
to (|I3l) . Throwing off time-like part of the action we can go again from the set 
of 3-dimensional Wegner models to the set of standard 3-dimensional Ising 
models via the duality transformation. And after that we may come back to 
variables of the original theory. 

K 'a H * K 'a s xS% +n 

Pa xn 

th< = i a = (18) 

1 + la 



7c. 1-7? 



c 



7(7 -i i 1 7r 1 i fr 7 

1 + 7c 1 + la 

Here tilde & (prime k') means duality transformations in three (four) dimen- 
sions, respectively. The original theory undergoes the phase transition at 
7 a ~ 1 and j T = j c . 

Similarly, let's consider the limit 7^ ~ (it corresponds 7 r ~ 1 for the 
original theory) as we have done already for the limit 7^ ~ 0. We can find 

1 1/N T 

the phase transition at 7 r ~ 1; 7<j = — xjnt- We may depict all these critical 

points in the plane 7 r <g) ^ a (fig-3). 

To avoid arising infinite constant in the partition function at 7^ = 1 and 
at 7 T = 1, we took 7^ and 7 T close to 1 but not equal 1 exactly. 
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ANOTHER APPROACH TO THE CONSIDERATION OF LIMITS: 

k t — > oo AND K a — > oo. 

There is also another method to find critical points 3 and 4 (fig. 3), different 
from the above duality transformation approach. Let us build up the effective 
action for the case K a — > oo. It is obvious that space-like plaquette variables 
in this case D nm = □<,. = 1. The gauge is completely fixed if we have a 
maximal tree, i.e. tree which will have a closed loop after adding one more 
link to it. 

We may choose the maximal tree in two dimensions like a snail as shown 
at the fig. 4. By analogy we may obtain such a maximal tree for 3 dimensions. 
That choice of the maximal tree provides that the equation n a = 1 will have 
only single solution: cr™ = 1. Really, if we consider the first plaquette, which 
have three link variables equal 1 due to the gauge condition and provided 
with D a = 1, the fourth link is to be equal 1 also. Now we can apply these 
considerations to the second plaquette in the snail and so on, getting finally 
cr" = l. As a result, only temporal links will survive in the action. 

- 5 = 7 T ^^cr (f,r)cr (f + n,r) (19) 

T xn 

We've got the set of 3-dimensional independent Isings with critical point 
at 7cr ~ 1; 7 r = 7 C for our original system. 

In the limit j T ~ 1 the plaquette variables n Qn = D T = 1. In the Hamil- 
tonian gauge (do = 1) cr n (t) = a n (t + 1) = a n are time independent. We 
come to the static 3-dimensional Wegner system 

- S = N T K a Vn{x)& m {x + n)a n (x + n + m)a m {x + m) (20) 

xn 

which transforms into the Ising model under duality transformations. 

1 + 7 e W 

So, this alternative method confirms the results obtained previously. 
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THE POTENTIAL BETWEEN PROBE SOURCES. 

We would like to estimate the connected correlation function (a x a x+R ) for 
the Ising model with different couplings in each direction within the spherical 
model. The crucial point is the following condition: 

^X>* = 1 (22) 



Then 

Z = J2 f + °° — e^ 3 -^^^*"^^ (23) 
^Jc-ioo 2iri 

where n' n = K m k the constant c is chosen to ensure the legitimacy of inter- 
changing the integration and summation order. It means that c is a line to 
the right of all ^-singularities. 

We can rewrite the partition function as: 



da 

z = ] 

where 



Z = j — e aN « ]T e-KAc-sXV (24) 



A x - X > = abt - £ n'J x x +n = J (a - £ < cos ^e^-^d 3 ^ (25) 

n=l n=l 



The correlation function ((J x (J x+ r) can be calculated as the derivative of gen- 
eration function over sources: 

(a x a x+R ) = (a a R ) = I ' da E e ^-^ x _ x ,* x , +Vx * x (26) 
and after shifting integration's variables we have 

(Wx+r) = 4r^T I r dae^- x ^> (27) 

g -i(f)R 



dqo dq R 

= A R 1 = I \— -d^ 



«o - Tr n= i < cos <f) n 
ao is the sadle point which is determined by the condition 

/ ^ T- = <"o> = 1 (28) 

J ao - E„=i «„cos0„ 
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At R n — > oo (j) n — > and the pole will be determined by 



2 n=l 

Introducing "symmetrical" variables 

c 



I 3 3 
o £ K n0n = -«0 + E < 



n=l 



0n = 



C = i 



2 («o - E O 

n=l 



we obtain 



/ \ ' t—<ri=\ /~T" 



(29) 



(30) 



(31) 



where r = i? n a n 

So, the potential between two probe sources depends on the choice of n. 



THE PHASE STRUCTURE AND LIMIT a Tj(T -> 0. 

We will make suggestions about the phase structure in the whole area of 
coupling constants and clarify the nature of the phases previously obtained. 
In the case K-r ~ the 4-dimensional system transforms into the set of in- 
dependent 3-dimensional subsystems with t = tj as mentioned already. The 
probe sources (the potential between them was calculated on dual lattice) 
correspond to the magnetic charges placed inside cubes of the original lat- 
tice. 

The production over space-like cube's plaquettes can be associated with 
magnetic field flux through the cube's surface 

n 

cube 

Bk 

and is not equal to zero when the probe source is placed in corresponding 
dual site. In the other words, the probe source of "electric" charge in the site 
of the dual lattice corresponds to the monopol ("magnetic" charge) of the 
original one. 



exp {const E B ■ ft} 

cube 

n^kmn^" mn (32) 
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As known, if for Wilson's loop C* in the plane [tx] of 4-dimensional dual 

lattice we'd fix the time t = to then the loop will pierce the plane [zy] in two 

points (monopol - antimonopol). Dirak's string which ties them together lies 

in the slice t = to []TD| . 

If the potential between probe sources in each slice will increase linearly 

i 

with R (in the region of coupling 7 T < 7^: = 7<f r ) then the average value 
of Wilson's loop (W) = Hf=o(s Xo s XR ) (when the slices are independent) will 
decrease exponentially according to area law. Average value of the corre- 
sponding t'Hooft's loop (t'H) must behave in the same way in the region 

la ^ la — 1+7c - 

T T 

(t'H) = H (m Xa , m XR ) = (W) du ai = II (s Xq s Xr ) ~ e~ aTR (33) 
t=l t=l 

It is obvious that parameters area (j T and j a ) falls to four sectors depend- 
ing on the behaviour of the average values of Wilson's and t'Hooft's loops. 



I 


7r > 7r 


lo < la 


(W) ~ e- aLc 


(t'H) - e~ a ' L o' 


II 


7r > 7r 


1° > la 


(W) - er aLc 


(t'H) - e- x '^c> 


III 


7r < 7r 


1° < la 


(W) - e- x ^ c 


(t'H) ~ e~ a ' L c 


IV 


It < It 


lo > la 


(W) - e~ ASc 


(t'H) - e- A '^' 



This picture covers all four types of possible behaviour of the averages 
under consideration which were found by G.t'Hooft jTTj from the commu- 
tation relations analysis. It seems impossible to "see" all four phases on a 
lattice with fixed asymmetry (k t = constK a ) including the symmetrical one. 
Our results are in good agreement with |TT| in the areas they studied on the 
symmetrical lattice (line 7 T = j a ) 



Up to now we treated the couplings 7 CT (7r) freely enough considering them 
independent. However, it is obvious that underlying constants ga{g T ) should 
depend on the lattice spacings through the renormalization group relations. 
This connection may make some areas of the square (0 < i a < 1; < 7 T < 
1) inaccessible. To find exact borders of the accessible area we should build 
the renormalization group relations on an asymmetrical lattice and put g a {gl) 
into 7o-(7 T ). In "naive" limit g a ~ g\. Introducing £ = ^ we get ^ = £ 2 . 
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Investigating the phase structure of our theory we were dealing with ef- 
fective couplings ft TjCr . Now we are interested in clarifying whether critical 
values fcj- (j are within accessible area of temperatures (3 1 = (a T N T ) 1 (i.e., 
(3 7^ and (3 ^ oo) or not. We should note that in the limit (a a , T — > 0) 
g 2 ~ g 2 ~ g 2 << 1 0. So, K T K a = ^p- >> 1 (shaded region at fig. 6). Say, 
at K(j — > 0; /'w ► oo as -r— and £ = — as -r— . This narrows the accessible 
parameters area. Moreover, at N T — > oo the points 1 and 4 move to point 
(7r = 1; 7a = 0) (see fig.6). 

Taking into account that in this of parameters 

i T ~ exp{-2iV r e~ 2 ^} ~ exp{-2/?- — -} (34) 

we may see the critical point ^ is accessible at finite temperature (3~ l = 
(a r A^ r ) _1 when 

#r = j^-T- (35) 



a T A T 

where A r = 



As /v,t depend not only on a T but also on a a (directly and via g a{r too), 
finding functional relations between K a<r and j3 seems impossible if we don't 



fix the connection between and A^ r 



2Na a N T _2N [3 
9r P Si a<jN 7 



K r = 2 - T T K<t = _ 9 I Ar (36) 



In the "naive" limit g 2 ~ ^ and 

Av r (a a N T )'' 



(37) 



The investigation of the SU(2) gluodynamics phase structure on asymmetri- 
cal lattice in the plane (3 x g 2 is beyond the frames of this paper and will be 
done later. 

The parameter £ = usually is chosen arbitrarily ( ^Hamiit = oo and 
^Euci = 1 )• So, if this parameter is not restricted with additional condition 
then changing it arbitrarily we may reach any point of the curve K a K T = const 
at any small g 2 , thereby crossing at least one line between phases (77 and 
IV). This says that thermodynamical quantities depend on £ and, moreover, 
the jump on this parameter is possible for some of them. It is commonly 
believed that changing the parameter £ should not result in any observable 
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effects. The renormalization group relations make possible excluding the 
dependence of observed quantities on £ in "naive" limit, thereby making 
lattice regularization with different £ equivalent [2], [T2], [T3j]. To save the 
independence of observed quantities from £ in "naive" limit, we chose 
precisely the same way as in [g, p~2|, |T3|1. In the approximation SU (N) ~ Z(N) 
we failed to create the renormalization group relations with A parameter 
being a function of £ which remove the lattice asymmetry parameter from the 
thermodynamical quantities keeping at the same time "naive" limit results. 
So, we suggest that lattice gauge theories need some additional condition 
which fixes £. 

Although all calculations were carried out for the Z{2) group, they can be 
fulfilled for the Z(S) group also, with small changes. There are reasons to 
hope the results for SU(2) gauge group will be similar to those for the Z(2) 
gauge group, at least within approximations of [|5j]. 

The authors are indebted for fruitful discussions with Prof. Adriano Di 
Giacomo. They are grateful to O.A.Borisenko for the critical notes. 
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Figure 2. 
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Figure 3. 

Points 1-4 and 2-3 - dual symmetrical by pairs. 
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Figure 4. 
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Figure 5. 

/ - deconfinement of electric and magnetic charges 
77 - magnetic confinement, electric deconfinement 

III - electric confinement, magnetic deconfinement 

IV - electric confinement, magnetic confinement 

V - we cannot investigate this region analytically 
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